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Abstract: We present a study of the effect of an electric field on the binding 
energy of a shallow hydrogenic impurity in GaAs / Ga>i_ x Al x As quantum well 
wires. The wire is considered to be of length L and radius R and the electric 
field F is applied along the z— axis, the axis of the wire. 



1 Introduction 

The development of the epitaxial crystal growth techniques such as molec- 
ular beam epitaxy and metal-organic chemical vapor deposition has made 
the growth of the quasi-two-dimensional (quantum well) , quasi-one- dimen- 
sional (quantum well wire), or quasi-zero-dimensional (quantum dot) become 
possible [1,2,3]. 

Several studies have investigated the optical and transport properties 
based on the calculated electronic structure of these systems in the presence 
of shallow impurities (see for example references [4-9]). 

Regarding the quantum well wire(QWW), the binding energy of hydro- 
genic impurities in a QWW has been calculated as a function of the wire 
radius in case of no fields applied[7,8] . Also, Branis et al [10] considered the 
case of a magnetic field applied parallel to the wire axis. They found that 
for a given value of the magnetic field, the binding energy is larger than that 
of the zero-field case. 

In this work we have considered the case of an electric field applied along 
the wire axis, we have calculated the binding energy of the ground state of 
a hydrogenic impurity located at the wire axis. A variational approach has 
been utilized with a trail wave function containing a hydrogenic part. The 



other part of the variational function which represents the wave function 
in the absence of the impurity has been calculated analytically in terms of 
Bessel and Airy functions. 



2 Electron Eigenstates 

Consider a quantum well wire of length L and radius R (of dimension compa- 
rable to de Broglie's wavelength), with an electric field of strength F applied 
along the axis of the wire . The Hamiltonian for such a system is given by: 

H=^-V 2 + \e\Fz + V(p,z), (1) 

where 

{0, p< R and \z\ < \ 
(2) 
oo, p > R or \z\ > |f, 

We have used the cylindrical coordinates such that the origin is taken at the 
mid point of the wire axis, which is considered as z-axis 

The Schrodinger equation takes the form 

1 d , dip, 1 d 2 ip d 2 ip 2m e \e\F —2m e E 

~pd~ P [p W + 7W ^~ z ^ = fr 2 ^ {) 

Assuming that 

i> = $(0)C(^(p), (4) 

and noting that the problem is rotational invariant ($(0) = e tm< ^), then Eq. 
(3) takes the form 

1 d , cfiftx m 2 ld 2 ( 2m e \e\F _ -2m e E 

W P Tp^ p ltp~> ~~ P Y + ~cjh? h 2 z ~ k 2 ' ( ' 

The radial part is 

9*(p) = J m (Kmp), (6) 
where J m is the ordinary Bessel function . 

Accordingly the axial part ((z) is found to satisfy the equation 

ld 2 ( 2m e \e\F 2m e E 2 

z + — k nm = 0. (7) 



(dz 2 h 2 h 2 
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To solve this equation, we assume that 

* = G*(G*-*2^! + »£,.), (8) 

where 

G (9, 
Substituting from Eq. (8) into Eq. (7) we obtain the equation 

^ = iC, do) 

which is Airy's equation The solution of Eq. (10) is a linear combination of 
the Airy functions Ai(z) and Bi(z). (Cetina and Montenegro [12])Thus, the 
final form of if)(r) is taken as 

iinUD = Ne m tj m (k nm p) \Bi{z) - !tt&Ai{z)\ , (11) 

where £j = z (z = The form of £(z) in Eq. (11) has been chosen so 

that ip vanishes at z — = ^ L . We still need to satisfy the boundary conditions 
(Vanishing of the wave function at the boundaries) 

V(p = R)=0, i/>(z = |) = 0. (12) 

The first gives 

Jm(k nm R) = 0. (13) 

For the ground state (m = 0, n — 1), which is considered in this work, the 
above equation implies 

k w R = 2.4048, (14) 

where the R.H.S. is the first root of Jo(x) = 0. Also, the second condition in 
Eq.(12) yields 

Bi^ h )Ai^ t ) = BH&Aifa), (15) 

where ^ = z (z = |f). 

The energy states can consequently be obtained by solving the transcen- 
dental Eq. (15). The energy of the ground state will be denoted by E . 
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2.1 Impurity Eigenstates 

We consider the problem of a hydrogenic impurity located at p = z = 0. The 
trial wave function for the ground state is taken as 

ij = NJ (k w R)(Bi(z) - ^M Ai{ ~ z))e -^V^. (16) 



The normalization condition leads to 
N~ 2 = 2tt 



R ,# 



p=0 J z 



-L 



dpdz P J 2 (k 10 R)(Bt(~z) - |||2^(S))V 2 V^ 

(17) 



The potential energy (V) is 



-e 2 



(V) = ^|__=|^> 



eo\/p 2 + z 2 ' 



V e > J p= J z== L r y^T^ 2 
x ^Bi(z)-^^Ai(z)y e - 2X ^^. (18) 

Putting the potential energy in unit of Ryrderg R b and length in unit of Bohr 
radius ag, (V) takes the form 



,/p=0 Jz==± 



x j 2 (w)(^(5)-§m^(5)) 2 . (19) 

To calculate the kinetic energy (T) we first calculate V 2 -0 as 

V> = 01 02, 

where 

^ = e -Wp > + z2 

and 

02 = J (k 1Q p)(Bi(z) - *^Ai{z)). 
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Hence 

V 2 ^ = (V 2 Vi)^2 + ^i(V 2 ^ 2 ) +2(Wi)(V^ 2 ), (20) 

but 



(V 2 - Gz)fr = -^% 2 (21) 



and 

2A 



V 2 ^=(A 2 --^=W (22) 
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We then evaluate the term 



2(^> | W>iW> 2 ) = d(t)dzdpp{V '^1)^2)^2* 

= \\ d<$>dzdp p^ 1 ^ 1 ■ ^ ? ^ 2 



<9p <9p cfe cfe y 

+ / # ^ PP (_^ f , f )M (23) 

Integrating the first term by parts with respect to p and the second term by 
parts with respect to z and using the boundary conditions (vanishing of ^ 2 
at the boundary) we get 

2{i/j I V^iW> 2 ) = -2A 2 / d(f)dzdp p^ 2 + 2A / d^dzdp—jJ==^\i)\. 

J J V P 2 + -2 2 

(24) 

Using Eq. (21),Eq.(22), and Eq. (24), we get 



z£<*|V» = * + (25) 

Putting Eq. (25) in units of Rydberg R b , and length in units of Bohr 
radius a# and adding to (V), the total energy is 



A A 

(H(R,L,F)) = E + \ 2 + 1I (26) 



where 



/ / ^ ^ = J 2 (^ioP)fa(^)- 7 ^(^(^)) • (27) 



Minimizing (H(R,F,L)) with respect to A, we obtain the total energy, and 
the binding energy is thus 

a A 

E b = E - (H(R, F, L)) = -A 2 - ^ (28) 

d\ 

3 Results and Discussions 

The variation of the binding energy with the electric field strength F and 
the wire radius R are given in figures (1) and (2) respectively. The results 
displayed in Figure(l) show that for a wire of fixed dimensions L and R 
the binding energy decreases as the strength of the applied electric field 
increases. This is due to the fact that the increase of the electric field causes 
the electron to be less confined to the impurity and accordingly reduces the 
binding energy. Also the results in figures (1), (2) indicate that for a fixed 
electric field the binding energy increases as the wire radius decreases until 
it diverges as R — > 0. The divergence of the binding energy occurs since the 
infinite confining potential of the well forces the electron and the impurity to 
become very close as R — > 0. The same behavior has been previously pointed 
out in the absence of the electric field (ref[7]). 

4 Conclusions 

The binding energy of the hydrogenic donor impurity decreases as the strength 
of the applied electric field increases. It behaves with the variation of the 
wire radius in a similar manner to the case of no electric field applied. 
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Figure 1: Variation of the binding energy in Rydberg (i?^) units as a function 
of the electric field strength F (kV/cm). The wire length L = 3 a^,while three 
values of the wire radius are considered 0.50 ag, 0.75 ag, and 1.00 a B . 
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Figure 2: Variation of the binding energy E b in Rydberg i? fe units as a function 
of the wire radius R in Bohr radius a^units, where, F = lOkV/cm, and 
L = 3ob. 
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